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Inspired by the ubiquity of composite filamentous networks in nature we investigate models of 
biopolymer networks that consist of interconnected floppy and stiff filaments. Numerical simulations 
carried out in three dimensions allow us to explore the microscopic partitioning of stresses and strains 
between the stiff and floppy fractions c s and Cf, and reveal a non-trivial relationship between the 
mechanical behavior and the relative fraction of stiff polymer: when there are few stiff polymers, 
non-percolated stiff "inclusions" are protected from large deformations by an encompassing floppy 
matrix, while at higher fractions of stiff material the stiff network is independently percolated and 
dominates the mechanical response. 
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The basic design of most structural biological materials 
is that of a crosslinkcd mcshwork of semiflexible protein 
polymers. The mechanical properties of these biomateri- 
als are biologically highly significant [l], Q ■ Understand- 
ing these properties at the bulk or continuous level is not 
sufficient: biological entities like cells, motor- and sensing 
proteins experience, manipulate and interact with these 
polymer networks at single-filament lcngthscales and are 
therefore intimately aware of the discrete nature of these 
materials. Another complication arises when considering 
that most structural biomatcrials are in fact composites: 
bi- or polydispcrsc mixtures of different protein polymers. 
The extracellular matrix (ECM) consists of a mixture of 
stiff collagen and flexible elastin filament (bundles), and 
the relative abundance of these two greatly affects me- 
chanical properties Q . A more specific example that de- 
rives much of its biological function from the side-by-side 
deployment of mechanically vastly different filaments is 
articular cartilage - a complex, partially orded compos- 
ite containing typc-II collagen and proteoglycans as its 
main structural components [4| ■ Composite physics may 
be at play in single-component networks: coexistent and 
interlinked single fibers and fiber bundles determine the 
mechanical properties of actin gels and actin-filamin net- 
works [H . The interplay between stiff and floppy elements 
goes far beyond simple property mixing: The network 
of relatively floppy f-actin and intermediate filaments is 
believed to be nonlinearly stiffened by the rigid micro- 
tubules, and experiments have hinted at significant ten- 
sional forces in the cellular actin The cell cytoskele- 

ton, that is built up from microtubules, actin filaments 
and intermediate filaments, is yet another striking exam- 
ple of a composite network. 

Significant effort has been devoted to model systems 



of homogeneous and isotropic single- component networks 
of biopolymcrs, such as f-actin and collagen 8 16]. 
The single filaments that constitute these networks can 
be described by the semiflexible wormlike chain force- 
extension curve, where extension requires that thermal 
fluctuations of the filaments be suppressed leading to a 
steep and non-linear increase in the force. Compression 
requires considerably smaller forces that become constant 
in the Euler buckling limit 0, [i| . Networks of such fil- 
aments show highly non-linear strain-stiffening and neg- 
ative normal forces under shear 111 ]. Recent theo- 
retical studies and simulations have also underlined the 
importance of non-affine bending deformations in these 
networks 12, 14-3]. Models applying a similar method 



to composite biomaterials have only recently begun to 
emerge [It], EH and have focused on bulk behavior. 

In this Letter, we report the results of a series of numer- 
ical experiments of two-component networks of biopoly- 
mers to determine the relationship between composition 
and mechanical properties, both on the single-filament as 
well as on the bulk level. Furthermore, we compare our 
results to a theoretical model. 

Our networks consist of long filaments that are per- 
manently crosslinkcd. These crosslinks force a binary 
bond between two filaments, without angular preferences. 
The filaments are described by the semiflexible worm- 
like chain model [H, Q. The segments (parts of fila- 
ments in between two crosslinks) have degrees of free- 
dom in addition to the curvilinear length and the lo- 
cations of both ends, that are integrated to obtain an 
effective Hamiltonian. This effective Hamiltonian gives 
the energy for each network configuration, characterized 
by the network topology, the positions of the crosslinks, 
and the link-to-link separation of the filaments between 
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FIG. 1: Macroscopic properties of the networks as a function of the fraction c s of stiff filaments in a network, (a) Shear modulus 
K as a function of shear 7, normalized by the initial shear modulus Koj of single-component networks of floppy filaments. 
The different curves represent the stiffness of networks with c 3 = 0.0,0.56, 1.0 (from bottom to top), at a fixed persistence 
length ratio R p = 16. Some data points (< 1%) lie well outside the curve; these are indicated by the symbols. These outliers 
occur due to local reorientations, cf. (b) The normalized initial shear modulus as a function of c s , for networks with 

R p = 64 and 16. For comparison, we also plot curves corresponding to a linear scaling of the shear modulus with c s , given by 
Ko(c s ) = CsKo,s + (1 — c s )Ko,f, and a linear scaling of the compliance with c s , given by l/Ko(c s ) = c a /Ko, a + (1 — c s )/Kqj. 
(c) Critical shear 7 C , defined as the shear at which the shear modulus is twice the initial shear modulus, as a function of c 3 . 
The curve is drawn as a guide to the eye. 



crosslinks. Starting from a random, isotropic network 
consisting of crosslinks and segments, we apply a large 
number of Monte Carlo moves which alter the network 
topology such that filaments with a persistent direction- 
ality along segments are formed. At this point, we des- 
ignate filaments to be either stiff or floppy by assign- 
ing to each segment a persistence length and an equilib- 
rium backbone length. We then further relax the con- 
figuration by applying new Monte Carlo moves. All our 
networks have periodic boundary conditions and contain 
1,000 crosslinks. Their lateral sizes are determined by 
the condition of zero pressure. A detailed description 
of this approach is presented in 14|. Our networks are 
characterized by the following set of parameters: the per- 
sistence length £ p of the stiff filaments, the stiffness ratio 
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the average filament length L, the av- 
erage distance between crosslinks along a polymer's back- 
bone £ c and finally the relative fraction of stiff filaments 
c s . In the present work, we examine the c s -dependence 
of the mechanical behavior. We restrict ourselves to a 
biologically relevant region of parameter space: the per- 
sistence length £ Pi fl op py of the floppy filaments and the 
crosslink distance £ c are of comparable magnitude. On 
average, each filament is crosslinked six times (L = 6l c ). 
The ratio of the persistence lengths of stiff and floppy 
filaments R p is chosen to be 16 or 64. While this ratio is 
smaller than that for collagen/elastin (R p w 100) 0] or 
microtubules/f-actin (R p > 200) [3, it is large enough 
to capture the qualitative behavior of such composite net- 
works. Unless otherwise stated, all data shown represent 
the averages of nine network realizations. 

Our key findings are summarized in Fig. [TJ In Fig. [TJi 
we plot a 50/50 composite: rather than averaged, the 
mechanical behavior is bimodal - approaching the fully 



floppy system at low strains but, at finite strains, resem- 
bling the fully stiff network. We stress that this type of 
response can only be achieved in a composite. Even the 
linear behavior does not interpolate simply between stiff 
and floppy: Fig. [T]d shows that at low to intermediate 
Cg, the modulus is quite insensitive to c s , but rises very 
quickly at higher c s . Fig. [TJ;, finally, reinforces the point 
of Fig. [TJi: although the effects of adding stiff polymer 
are hardly noticeable in the linear elastic behavior, their 
effect on the nonlinear behavior is felt much earlier. The 
critical strain (7 C ) for the onset of the nonlinear regime 
reacts immediately to the addition of stiff material, but 
saturates at a point roughly coincident with the rise of 
the linear modulus. 

The qualitative picture that emerges at small c s is 
one of a floppy matrix encompassing isolated stiff fila- 
ments, or non-percolated clusters. Intuitively, the initial 
insensitivity of the linear modulus Kq to the addition of 
stiff material makes sense: deforming the stiff filaments 
requires higher energies than deforming the softer ele- 
ments, and therefore the low-energy modes of the system 
favor straining the floppy elements over the stiff ones. 
As long as stiff filaments do not form an independently 
load-bearing subnetwork, these low-energy modes exist 
and are compatible with the bulk deformation. 

This interpretation is confirmed by an examination of 
the microscopic deformation field, characterized by the 
non-affinity parameter A = (|x a fj — x| 2 )/7 2 (Fig. [2Ji). 
This parameter quantifies the deviation of the local de- 
formations, x, from a homogeneous affine deformation 
field, x a fj. As shown in previous numerical work (lH | 
and in experiments pol ]. the non-affinity generally in- 
creases with increased stiffness of the filaments, as bend- 
ing deformations are more important for the network 
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FIG. 2: (a) The non-affinity at zero shear, divided by the initial non-affinity of a network with c s = 0. The curves are drawn 
as guides to the eye. (b) Average forces in the floppy and stiff filaments during deformation, shown by the solid and dotted 
curves, respectively. The black curves represent the average force of a network with c s = 0.56 during shear. For comparison, 
we plot the average forces in single-component networks, c a = 1.0 and c s = 0.0 (grey curves). As indicated by the arrows, the 
curves for the average forces in the composite networks are shifted along 7 with respect to the curves for the single-component 
network, (c) The scaled initial stiffness as a function of c a , obtained by the floppy-mode model for R p = 16, 64, 1000, 00 (from 
top to bottom). For comparison, the data from simulations are given by the symbols. 



response of stiffer filaments. Indeed we find that the 
non-affinity is minimal for purely floppy networks and 
rises roughly linearly with the addition of stiff material. 
Such a linear increase represents the generic behavior of 
low-density (stiff) inclusions that independently perturb 
the deformation field of their surrounding (floppy) ma- 
trix plj . These additional non-affine deformations bring 
the floppy filaments closer to the nonlinear part of their 
force-extension relation, giving rise to the decrease of the 
critical strain j c as displayed in Fig. [lb. 

In the nonlinear regime the inherent stiffening of a sin- 
gle semiflcxiblc polymer makes the distinction between 
floppy and stiff fractions highly strain-dependent, with 
the ratio of their nonlinear moduli tending to unity in 
the high strain limit. This suggests a self- matching be- 
havior at finite strain: the floppy network stiffens up to 
the point where its modulus matches that of the stiff net- 
work. Beyond this point, the entire meshwork behaves 
as a nearly monodisperse system of stiff filaments. This 
effect is the origin of the behavior in Fig. [TJt: at high 
strains, the entire system is ultimately forced to couple 
to the stiffer deformation modes. 

This mechanism of stiffness matching is illustrated in 
Fig. r2b> which shows the average forces in the stiff and 
floppy filaments during deformation. By comparing with 
the one-component networks (grey lines) we can define a 
strain shift ^7: For given network strain 7, the filaments 
in the composite behave as if they were strained up to 
7 + ^7. Apparently, the effective strain on the floppy 
filaments is much larger than that on the stiff filaments. 
Equivalently, high forces in stiff filaments are suppressed, 
at the cost of increased forces in floppy filaments. 

Interestingly, this load-partitioning persists even at 
zero strain, where stiff filaments arc, on average, com- 
pressed while floppy filaments are stretched out. This 



stretched/compressed ground state is tantalizingly rem- 
iniscent of tenscgrity states Q- Apparently, dense 
crosslinking restricts relaxation of the network, and the 
absolute minimum of mechanical energy cannot be at- 
tained. There is, therefore, always a finite amount of 
residual elastic energy. This suggests that such force dis- 
tributions may not be a deliberate design principle but 
rather are the necessary byproduct of polydispersity in 
filamentous composites. 

The picture of a floppy matrix embedding stiffer in- 
clusions breaks down when the stiff filaments become 
independently rigidity percolated: the point where de- 
formation of the stiff elements becomes inevitable. We 
may estimate the percolation threshold by a counting ar- 
gument [22l | . Equating the number of degrees of freedom 
of the stiff filaments to the number of constraints due to 
crosslinks between stiff filaments gives (for L/t c = 6) a 
threshold c s = 0.56. This marks the transition from the 
low to the high c s regime and coincides roughly with the 
rise in the linear modulus Kq (Fig. [Tb>). Two separate 
obervations confirm the onset of stiff dominance: Firstly, 
c s = 0.56 is the point at which the non-affinity, which 
we attribute to the floppy matrix attempting to work 
around the stiff fraction, begins to plateau at the level of 
the bending dominated response of a purely stiff network. 
Secondly, the critical strain 7 C levels off around this same 
value of c s . In the range of stiffness ratios (Rp) accessible 
to the simulations the percolation is rather "soft" and 
represents a smooth cross-over phenomenon. The ap- 
proach towards the singular percolation limit, R p = 00, 
has for example been studied in simulations of mixed 
random resistor networks [23j . To address the analogous 



problem we compare our results with theoretical consid- 
erations, in which the parameter R p can be tuned to ar- 
bitrarily large values. The "floppy-mode" theory [l2| has 
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recently been shown to capture quite well the elasticity 
in one-component isotropic 2J, |25[ as well as anisotropic 
networks [26j . Within this theoretical framework the cal- 
culation of the network elastic modulus is reduced to the 
description of a "test" filament in an array of pinning 
sites. The coupling strength to these sites, k, repre- 
sents the elastic modulus of the network and has to be 
calculated self-consistently. To generalize this model to 
the case of composite networks we use two different test 
chains with coupling parameters kf and fc,, representing 



floppy and stiff filaments, respectively [27|. The use of 
two different coupling strengths quite naturally takes into 
account the load partitioning encountered in the simula- 
tions. The network modulus, k = c s k s + (1 — c s )kf, is 
obtained by solving the two equations 
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where k ai — k s , kf with probability c s and 1 — c s , respec- 
tively. The two energy contributions on the rhs of Eq. ([1]) 
reflect the competition between the bending energy of 
the (floppy or stiff) filament, W^ s , and the energy due 
to deformation of the surrounding medium by displacing 
the pinning sites (located at arclength position Si along 
the filament). The nonlinear entropic stretching elastic- 
ity is not included in these equations. The minimization 
is to be performed over the contour of the filament, y(s), 
the angular brackets specify the disorder average over the 
network structure. 

Fig. [2b displays the results from this calculation for 
various stiffness ratios R p , showing a sharp percolation 
transition in the limit R p — ¥ oo. The model compares 
well with the simulation data, even though the stretch- 
ing elasticity is not accounted for. This indicates that 
the bending stiffness is likely the dominant factor in de- 
termining the rise of the linear elastic modulus, in agree- 
ment with the proposed mechanism of load-partitioning 
and the observed increase of the non-affinity 

In conclusion, our results demonstrate that the me- 
chanical behavior of filamentous composites is consider- 
ably richer than the simple proportional mixing of prop- 
erties. The fact that the floppy and stiff networks are 
physically linked causes a strongly nonlinear coupling 
between the strain fields which deeply affects compos- 
ite mechanics. This may explain the ubiquity of com- 
posites in structural biological applications: slight vari- 
ations in composition cause large changes in mechanical 
behavior. This high susceptibility makes the compos- 
ite architecture an attractive motif for biological regu- 
lation. Likewise, the "best of both worlds" aspect may 
be exploited by Nature: composites combine the initial 
softness of their most compliant components with the ul- 
timate toughness of the stiffest elements. This greatly 



enhances the stiffness range of nonlinearly clastic ma- 
terials. Moreover, composites do so in a manner that 
could never be attained in monodispcrse materials, since 
linear and nonlinear properties of composites are deter- 
mined by two physically different materials and therefore 
may be independently varied. This possibility of inde- 
pendently tuning the linear and nonlinear behavior also 
has considerable potential for the design of biomimetic 
or bio-inspired synthetic materials and deserves further 
exploration. 
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